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Comment on "Why Do Gallium Clusters Have 
a Higher Melting Point than the Bulk?" 

The computational work Q], motivated by recent ex- 
periments |31 on the heat capacity measurements of small 
Gaj^g and Ga4g clusters, claimed that the observed broad 
peak in the heat capacity represents melting of the small 
size clusters and made a strong point that due to the spe- 
cial character of the chemical bonds, these clusters, con- 
trary to all expectations, melt at temperatures higher 
than the corresponding bulk material. In order to un- 
derstand mysterious "higher-than-bulk melting temper- 
atures" in small galhum clusters, Gaiy and Gaia, they 
utilized the powerful machinery of the density functional 
theory (DFT) molecular dynamic (MD) simulations in 
a form of the isokinetic Born-Oppenheimer MD, using 
ultrasoft pseudopotentials within the LDA. The specific- 
heat curve, calculated by the multiple-histogram tech- 
nique, indeed showed the peak in the heat capacity to 
be well above the bulk melting point of 303 K, viz., 
around 650 and 1400 K for Gan and Gaia, respectively. 
The "higher-than-bulk melting temperatures" were at- 
tributed mainly to the covalent bonding in clusters, con- 
trasting the covalent-metallic bonding in the bulk. 

In our Comment we show that the peak in the heat 
capacity is in fact a well known generic behavior of fi- 
nite size systems usually referred to as Schottky anomaly. 
Thus the connection of the peaks calculated in 01 (and 
also in the earlier works to the real world is ques- 

tionable. 

The thermodynamic behavior of a finite system con- 
sisting of N particles is well discussed (as an exercise) 
in the R. Kubo textbook |^, Chapter 1, Example 4 on 
pages 38-41. When a system contains a substance having 
the excitation energy AE, the specific heat is given by 
the formula (9) of Example 4: 
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Shown in Figure 1 is the corresponding peak. The text- 
book exercise refers to the system of noninteracting parti- 
cles. To see qualitatively what do interactions contribute 
to Schottky anomaly, one can consider, for example, a 
finite system of interacting spins. To avoid unnecessary 
complications we take the simplest analytically solvable 
model for a spin chain closed into the ring, i.e. the one- 
dimensional Ising model with the periodic boundary con- 
ditions 0. The specific heat per spin 
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(J is the coupling constant) is shown in Figure 1. Note 
the shift of the position of the maximum to higher tem- 
peratures and growth of the maximum magnitude for the 
smaller cluster size exactly as presented in Figs. 4 and 5 
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FIG. 1: Left: The Schottky specific heat. Figure 1.9 from 
Right: Specific heat Cv per spin of the Ising spin chain with 
= 17 (sohd line) and A'' = 13 (dotted line) spins. 


of 1]. This kind of the size dependence is indeed typi- 
cal for the Schottky anomaly. Interactions thus can only 
slightly change its shape, but the peak itself remains in- 
tact. This concludes our demonstration that the broad 
maximum in the heat capacity observed in 2] should be 
attributed rather to Schottky anomaly, than to cluster 
melting. 
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